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We present predictions at NLO accuracy in QCD for top-quark pair production induced by an 
anomalous chromomagnetic dipole moment of the top quark. Our results are obtained for total as 
well as fully differential cross sections, including matching to parton shower simulations. This process 
is expected to provide the most stringent direct limits on top-quark chromomagnetic dipole moment. 
We find that NLO corrections increase the contribution from the dipole moment by about 50% at 
the LHC, and significantly reduce the renormalization and factorization scale dependence. Using 
the NLO prediction, we update the current limit from the Tevatron and the LHC measurements. 
Apart from total cross section, we also study other observables relevant for LHC phenomenology. 


I. INTRODUCTION 


The top quark is expected to play an important role 
in new physics searches, due to its large mass and strong 
coupling to the electroweak sector. Strategies to search 
for new physics effects in the top sector can be broadly 
divided into two categories. In the first category, we 
search for new resonant states, such as tt resonance and 
top partners. In the second category, new states are as¬ 
sumed to be too heavy to be directly produced, and their 
indirect effects are searched for in top-quark couplings. 

Currently, no new states have been discovered, and 
exclusion limits have been placed, up to around several 
TeV scale for many new particles in either complete or 
simplified models [T]. On the other hand, in the sec¬ 
ond case the interaction of the top quark is becoming an 
ideal probe to new physics. On the experimental side, 
the millions of top quarks already produced at the LHC 
together with the tens of millions expected in the com¬ 
ing years will move top physics to a precision era. Many 
detailed and accurate information on various top-quark 
properties have been collected, and more will come. On 
the theory side, accurate SM predictions are also avail¬ 
able, in general at next-to-next-to-leading order (NNLO) 
in QCD and next-to-leading order (NLO) in electroweak 
for inclusive observables and at NLO in QCD for more 
exclusive ones. All of these provide the opportunity to ex¬ 
tract or to constrain different anomalous top-quark cou¬ 
plings. In this context, theoretical predictions including 
QCD radiative corrections to anomalous top-quark cou¬ 
plings will be necessary for extracting precise and reliable 
limits, as leading order (LO) predictions in hadron col¬ 
liders are often not reliable and suffer from large scale 
uncertainties. This has motivated a significant activ¬ 
ity dedicated to providing the NLO QCD corrections to 
top-quark processes involving anomalous couplings, or 
higher-dimension operators EHH]. However, NLO pre¬ 
dictions involving anomalous top-quark interactions are 
still far from complete. 

In this work we focus on the chromomagnetic dipole 
moment (CMDM) of the top quark in tt production. The 
cross section of tt production is one of the most accurately 
measured observables in top physics, and the effect of an 


anomalous CMDM has been investigated in many studies 
[TBH2H- To the best of our knowledge, the contribution of 
top CMDM is known only at LO. The goal of this work is 
to promote it to NLO including fully differential produc¬ 
tions and matching to parton showers, and investigate 
its impact on the current limits of top-quark CMDM, as 
well as other observables, with either stable or decayed 
ti system. We shall mention that, apart from CMDM, 
the top quark can also have an anomalous chromoelec¬ 
tric dipole moment (CEDM), which in this work we will 
not discuss. The reasons will be explained in the next 
section. 


At first glance, one might expect deviations induced 
by anomalous top CMDM to be small, and therefore ra¬ 
diative corrections of these contributions to be a higher- 
order effect. However, in tt production at the LHC, the 
K factor from NLO QCD corrections is about 1.5 in the 
SM, and is numerically not a higher-order effect. If a 
similar K factor applies to the CMDM contribution, it 
will be important to know the NLO correction to the 
CMDM, so that a more accurate and stringent limit can 
be obtained. Moreover, this statement is too naive, since 
the issue is both on the accuracy, i.e. the central value, 
and on the precision, i.e. the uncertainties of a prediction. 
LO predictions for processes at Hadron colliders always 
suffer from large uncertainties due to scale variation, and 
the NLO predictions are expected to significantly reduce 
these uncertainties and thus provide a more reliable es¬ 
timation of the possible range of the anomalous CMDM. 
Finally, in many cases NLO corrections can have an im¬ 
pact on kinematic distributions. Knowing the accurate 
differential cross section from the CMDM is therefore im¬ 
portant in measurements where the shapes of the distri¬ 
butions are used. We will show such examples, where 
LO prediction for the distributions does not provide a 
reliable description. 


The paper is organized as follows. In section [TT] we 
briefly discuss the theoretical background of top CMDM. 
In section uni we describe the framework of our calcula¬ 
tion and how it is implemented. Our results for total 
cross sections and limits are presented in section IV In 
section[V]we show several examples of exclusive distribu¬ 
tions. Section El is our conclusion. 
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II. THEORETICAL BACKGROUND 


The top-quark chromomagnetic and chromoelectric 
dipole moments, CMDM and CEDM, can be parame¬ 
terized by adding an effective term to the top-gluon cou¬ 
pling: 

C ttg = g s t^T A tG A + la* v (dy + id a 7s) T A tG A v 

( 1 ) 

where g s is the strong coupling, and G A U is the gluon field 
strength tensor, dy and dA in the second term represent 
the CMDM and CEDM of the top quark respectively. 

The CMDM of the top quark can arise from various 
models of new physics. The Yukawa corrections to gti 
vertex in two Higgs doublet model (2HDM) was first con¬ 
sidered in Ref. [16], while the supersymmetric QCD and 
electroweak corrections have been studied in Ref. uzhis]- 
Explicit expressions for CMDM in 2HDM and in minimal 
supersymmetric standard model were given in Ref. |20| . 
The top CMDM also arises quite naturally in composite 
models and technicolor models m ■ For a more general 
discussion of top CMDM in new physics scenarios we 
refer to Ref. [22]. Finally, a top CMDM operator may 
be loop-induced by operator mixing effects, from other 
higher-dimensional operators generated at higher scales. 
An example can be found in Ref. EDI- 

At the LHC, CMDM is mainly constrained by tt pro¬ 
duction. Direct limits have been derived by previous 
studies [25IiTT)l . However, the contribution of top CMDM 
has been known only at LO accuracy. Our aim is to pro¬ 
vide the NLO prediction, as well as to study its impact 
on the total cross section and various distributions. 

To go beyond LO calculation, a theoretical framework 
based on the dimension-six Lagrangian of the SM is re¬ 
quired. This framework contains a complete set of opera¬ 
tors satisfying the symmetries of the SM, i.e. the Lorentz 
symmetry and the SU (3)c x SU(2) L xU (l)y gauge sym¬ 
metries. It provides an unambiguous prescription for op¬ 
erator renormalization, and thus allows for a complete 
and consistent treatment of the higher-order corrections 
to the operators. The Lagrangian including dimension- 
six operators can be written as 


£eft = £-sm + ^ 


C 

A 2 


+ h.c. 


( 2 ) 


where A is the scale of new physics. In this work we work 
up to order 0{ A -2 ), as going beyond this order would 
require complete knowledge of dimension-eight operators. 

The top-quark CMDM in this framework is represented 
by a dimension-six operator 

Ota = ytgs (Qa^T A t) 4 >G a v , (3) 

where Q is the left-handed top- and bottom-quark dou¬ 
blet, t the right-handed top, <j> the Higgs doublet, and y t 
the Yukawa coupling of the top quark. <j> = ia 2 (/). This 
operator, after the electroweak symmetry breaking, takes 






FIG. 1. Representative tree-level diagrams of tt production 
with an effective vertex form the operator OtG- Black dot 
represents effective vertex from OtG- 


the form of the second term in Eq. ([!]). The relation be¬ 
tween dy and the real part of the coefficient of O t G is 
given by 


iv = (4) 

A“ 

The operator OtG contributes to tt production at tree 
level by modifying the standard gtt vertex, as well as 
inducing a new ggtt vertex, as shown in Figure [l] The 
effects of this operator in top-quark processes at LO in 
QCD have been discussed in Refs. EHEH- 

On the other hand, dA, the CEDM, corresponds to the 
imaginary part of C t G- In this work, however, we are go¬ 
ing to focus only on the CMDM. This is because the anal¬ 
ysis of the CEDM at NLO follows a completely different 
approach. As we have mentioned above, in an approach 
based on the dimension-six Lagrangian, we can only work 
up to order 0{ A~ 2 ), and thus only the interference be¬ 
tween the CEDM and the SM amplitudes can be in¬ 
cluded. At this order the contribution vanishes in tt pro¬ 
cess because of the CP-odd nature of the CEDM, unless 
one incorporates the decay of the top quarks. As we shall 
see, our work is based on the MadGraph5^\MC@NLO 
framework [ 55] . where the spin correlation and the off- 
shellness of the top quark pairs are simulated by using 
the MadSpin package [31], which is based on LO evalua¬ 
tion of the complete matrix element including top decays. 
Therefore it is not a suitable framework for the NLO cor¬ 
rections to the CP-odd effects of the CEDM, and we will 
leave the NLO analysis of the CEDM to future works. 
Throughout the paper, we assume CtG to be real. 

When going to NLO in QCD, one needs to take into 
account the operator mixing effects between OtG and 
other dimension-six operators that could give a contri¬ 
bution to the same process at tree level. For tt pro¬ 
duction, these operators are Oq = g s f ABC G Av G^ p G c p ->*, 
0(j>G = gl {$<!>) G a „G a ^ , and several four-fermion op¬ 
erators [ST]. It turns out that in tt production, the mixing 
from Otc to these operators is not relevant. First of all, 
OtG does not mix into Oq and four-fermion operators 
[30] , because O t G is essentially a dimension-five operator 
if the Higgs field always takes the vacuum expectation 
value, which is always true at the order we are working 
at. Second, O t G does mix into O^g [35], but such effects 
correspond to a 0{y 2 ) correction to the LO process, and 
therefore of higher-order. Finally, the operators O^g an d 
Oq do mix into OtG !3E;j however it is consistent to as- 
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sume that they vanish at all scales, given that they are 
not renormalized by OtG- Therefore, as a first step, in 
this work we assume O t G is the only non-vanishing op¬ 
erator, and neglect other operators. Note, however, that 
for a fully consistent phenomenological study, the com¬ 
plete operator set must be included. A first example of 
the global approach was presented for the flavor-neutral 
interactions of the top quark [37]. The NLO predictions 
for other operators will be left to future works. 



FIG. 2. Contribution of OtG operator to top-quark and gluon 
wave functions. Black dot represents effective vertex from 
OtG- 


III. FRAMEWORK AND IMPLEMENTATION 


In an NLO calculation one has to choose a renormal¬ 
ization scheme. Our scheme is consistent with Ref. [S]. 
For the SM part, we adopt MS with five-flavor run¬ 
ning in a s with the top-quark subtracted at zero mo¬ 
mentum transfer [38j . The bottom quark mass is ne¬ 
glected. Masses and wave-functions are renormalized on 
shell. The dimension-six operator OtG then gives addi¬ 
tional contributions to top-quark and gluon fields renor¬ 
malization, as shown in Figure [2] We find 


SZ^ 

Sm t 

SZ^ 

where 


SZ. 
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\£uv 


5m t ,sM — C t G —“ ttw " D 


Aotsirit ( J_ 1 


SZ. 


(a) 

2 ,SM 


-C tG 


7tA 2 
2 u s m 2 t 1 

7tA 2 £ Ejjv 


\ + q 

\£UV 3 


D.^Td + e)^ 
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(5) 

( 6 ) 

(7) 

( 8 ) 


and /i is the renormalization scale. In addition, the strong 
coupling counterterm, Z 9b , also gets a dimension-six con¬ 
tribution: 

cv m 2 1 

SZ g , = SZ 9s! sm + C tG ^D e - , (9) 

Ttl\ £jjv 

that is to say the top-loop contribution with the operator 
OtG is also decoupled from the running of a s , in the same 
way as in the SM. Finally, for operator coefficient we use 
MS subtraction. The counterterm of C t G is 

sz Cta = ^r(i + £)(^r (io) 


parameter. Therefore it is equivalent to redefine Ot$ as 
yf {<jy<j> —v 2 /2) (Qt)(j), and shift the renormalization of the 
dimension-four component of O t ^, i.e. — m 2 /A 2 yt(Qt)cf), 
to Eq. |6|. This is more convenient since O then com¬ 
pletely drops out from the calculation. 

Our calculation is performed using the Mad- 
Graph5_aMC@NLO framework [33] . The operator 
OtG is implemented in the UFO format s32| by using 
the FeynRules package [40]. Helicity amplitude rou¬ 
tines are generated by ALOHA HU- The evaluation of 
the loop corrections requires two additional pieces, the 
UV counterterms and the rational R2 terms which are 
required by the OPP technique [42]. The UV countert¬ 
erms are computed according to Eqs. <®> -Pt, while the 
R2 terms are generated by the NLOCT package [13] • The 
calculation is then automatically performed by Mad- 
Graph5_aMC@NLO at NLO accuracy, and matched 
to parton shower via the MC@NLO formalism [HI . 

Several checks of the implementation have been done, 
including the gauge invariance of all virtual contribu¬ 
tions, and the pole cancellation when combining virtual 
and real contributions. In addition, we checked that all 
relevant UV and R2 terms are correctly implemented, 
by computing individual diagrams with MadLoop [15] 
and comparing with analytical results obtained by using 
FormCalc and LoopTools [46] . 

IV. TOTAL CROSS SECTION 

In this section we give the NLO total cross section 
from OtG-, and place limits on its size, using available 
measurements from the Tevatron and the LHC. 

As mentioned above, we work up to 0(A~ 2 ), which 
means we insert in each diagram at most one effective 
vertex from OtG- The total cross section then becomes a 
quadratic function of C t c/ A 2 , 


This will lead to the running of C t G- 

One remark on Eq. (j6j) is in order. Naively if MS' is ap¬ 
plied to the complete set of dimension-six operators, one 
would expect that the dimension-six “Yukawa” operator, 
(f>)(Qt)<p, will be renormalized by O t G and 
will be providing the UV pole in the mass counterterm 
in Eq. ([6]). The remaining finite term, however, still needs 
to be subtracted by introducing the mass counterterm. 
Operator O t( f, does not have a physical effect in this pro¬ 
cess, as it only shifts the top-quark mass which is an input 


a - ct sm + 

The /3i term represents the contribution from OtG at or¬ 
der C>(A -2 ). The quadratic /3 2 term, on the other hand, 
does not have a physical meaning without a complete cal¬ 
culation at 0{ A -4 ), and needs to be dropped. However 
the size of this term can be used to gauge the range in 
which the approach itself is valid, or in other words, the 
expansion in 1/A converges. 
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ft 

LO [pb TeV 2 ] 

NLO [pb TeV 2 ] 

K factor 

Tevatron 

1 £1+0-66 (+41%) 

1 - O1 -0.43 (-27%) 

i si n+ 0 073 (+ 4 - 05 %) 

l.OlU_ 0 19 7 (_ 10i88 %) 

1.12 

LHC8 

r n t+ 17.3 (+34%) 
ou *'-12.4 (-25%) 

70 £0+9.26 (+12.7%) 

/ Z.OZ — 10 53 (_ 14>5 %) 

1.43 

LHC13 

1 £1 £+48.0 (+29.7%) 
101.0 36 2 (-22.4%) 

9*30 ec+29.0 (+12.1%) 
zoy,J - 31 .8 (-13.3%) 

1.48 

LHC14 

1 Q1 Q+55.6 (+29.0%) 

1 c/1 .o_ 42 2 (_22.0%) 

900 r»+33.6 (+11.9%) 
ZOO.U_ 3 g g (_ 13>1% ) 

1.48 


ft 

LO [pb TeV 4 ] 

NLO [pb TeV 4 ] 

Tevatron 

0.156 

0.158 

LHC8 

8.94 

11.8 

LHC13 

30.0 

43.2 

LHC14 

35.7 

51.6 


TABLE I. Values of /3i and /?2 at LO and NLO precisions 
for the Tevatron, and for the LHC at 8, 13, and 14 TeV. 
The respective K factors for the central values of /3i are also 
shown. 


To extract ft^, we perform the calculation with CtG 
taking different values: 0, ±1,±2, and fit the resulting 
cross sections to Eq. (111. Each run is performed with 9 


combinations of where fiR is the renormaliza¬ 

tion scale and /ip the factorization scale, each can take 
values ft/2, /i and 2/i, with the central value /i = rrit- 
This allows us to extract the scale variation of /3 1 . In our 
calculation m t = 173.3 GeV, and we use the NNPDF 2.3 
set of the parton distribution functions S3■ The values 
of flip, both at LO and NLO, are given in Table |T] for 
Tevatron, LHC 8 TeV, LHC 13 TeV and LHC 14 TeV 
runs. A significant improvement in the scale dependence 
can be noticed. A large K factor is found at the LHC. 
The sizes of ft implies that the effective approach is valid 
given that Cft/A 2 < 1 TeV -2 . Our LO results agree 
with Ref. [58] once we take into account scale variation 
and note the opposite sign convention of dy- 

With these results we can set bounds on the size of 
O t G using total cross section measurements. We replace 
the osM in Eq. by the most precise SM predictions 
at NNLO+NNLL accuracy in QCD, which are <Tgj+ = 
7.148 ± 0.218 pb and Og^ 0 = 244.9 ± 9.7 pb respectively 
for Tevatron and for LHC at 8 TeV [55]. We sum the scale 
and PDF uncertainties in quadrature and symmetrize the 
error around a central value. The combined measurement 
at the Tevatron (LHC) is = 7.51 ±0.40 pb (fTexp C = 
240.6 ± 8.5 pb) m El], where we have corrected for the 
top-mass difference using the prescription given in these 
references. The value of C t G can be extracted using 


CtG 

A 2 


^exp ^SM 


( 12 ) 


together with its corresponding uncertainty, given by 

i 

_2 , _2 / \ 2 2 
fc exp e SM / Cft \ 

ft exp - CT SM ) 2 \ ft / 

(13) 

where the experimental error e exp and the theoretical 
NNLO error egM are summed in quadrature, and ep 1 is 


SC tG 

^exp ^SM 

A 2 

ft 



LO [TeV -2 ] 

NLO [TeV -2 ] 

Tevatron 

[-0.33, 0.75] 

[-0.32, 0.73] 

LHC8 

[-0.56, 0.41] 

[-0.42, 0.30] 

LHC14 

[-0.56, 0.61] 

[-0.39, 0.43] 


TABLE II. Limits on CtG I A 2 . The corresponding lim¬ 
its combining Tevatrion and LHC8, in terms of dy, is 
[-0.0099,0.0123] at LO and [-0.0096,0.0090] at NLO (note 
the opposite sign convention of dy in [58]). For LHC14 we 
assume a 5% experimental error. 


the (symmetrized) error of /?i due to scale variations. We 
assume no correlation between the SM NNLO prediction 
and ft. One could also add in Eq. (131 a term represent¬ 
ing the error from the missing 0( A --1 ) terms, which can 
be estimated using ft(Ctc;/A 2 ) 2 , but the changes in the 
limits are negligible. 

We show the 95% CL allowed region for C t G in Ta- 
ble[TT] The improvement of NLO calculation for Tevatron 
is mild due to the small K factor. In the LHC cases the 
allowed range is significantly reduced at NLO. We also 
give the expected limit at the LHC 14 TeV run, assuming 
an experimental error of ±5%. 


V. DISTRIBUTIONS 

Our calculation is implemented via the Mad- 
Graph5_aMC@NLO framework, therefore simulation 
of any observable is automatic. In this section we present 
a few representative distributions of variables of par¬ 
ticular relevance for LHC phenomenology. To simulate 
parton shower we have used the Herwig 6 code 1)21 . 
Other shower programs are also available, including Her- 
wig++ [53] and Pythia 8 [55] , 

A. Stable top quarks 

We first look at kinematic observables constructed 
from stable top-quark pairs. In Figure [3] and Figure [4] 
we show the invariant mass of the top anti-top system 
and the transverse momentum of the top-quark for the 
LHC at 8 TeV, in each case with the differential K factor 
displayed in the lower panel. The contribution from OtG 
is extracted by generating event samples with CtG = ±2 
separately and taking the difference, in order to get rid 
of the quadratic terms in C+. These observables can 
serve as discriminators in case any deviation form the 
SM is observed, and will be useful in determining the 
type of new physics [32]. One can see that the NLO 
computation reduces the scale variation. The differential 
K factor is not a constant and drops at higher scales, 
however in both distributions we observe that the K fac¬ 
tor of the OtG contribution is similar to that of the SM 
contribution, so using the SM K factor to rescale the LO 
event samples from OtG can be a good approximation 
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FIG. 3. Top quark pair invariant mass distribution at LHC 8 
TeV. 



FIG. 4. Top quark transverse momentum distribution at LHC 
8 TeV. 


of the complete NLO result. Note that we have chosen 
Ctc/A 2 = 1 TeV -2 for convenience, even though this 
value is already excluded by the current limits. One can 
always rescale the curve to get corresponding result for 
any other value of CtG, as we have already removed the 
quadratic dependence on CtG- 

In Figure [5] we show the top quark pair invariant mass 
distribution for LHC at 14 TeV, at high mass region 
above 1 TeV. In this calculation we have set y = 1 TeV. It 
has been suggested that kinematical cuts on the ti invari¬ 
ant mass can improve the sensitivity to access the Otc 
operator, both in the angular distributions m and in 
the total cross sections [38] , This is because higher mo¬ 
mentum transfer is favoured by the dipole structure of 
OtG- However, using NLO (or NNLO) prediction for the 



no cuts 

m t t > 1 TeV 

m t t > 2 TeV 

K (SM) 

1.49 

1.16 

0.77 

K {Ota) 

1.49 

1.14 

0.69 

O tG (LO)/SM(LO) 

0.32 

0.28 

0.29 

O tG (NLO)/SM(NLO) 

0.32 

0.28 

0.26 

OtG (LO)/SM(NLO) 

0.21 

0.24 

0.37 


TABLE III. K factor and signal excess u(Ot G )/cr(SM), with 
no cuts, and with cuts m t t > 1 TeV and 2 TeV. y = m-t- 
Both LO and NLO give similar results for the signal excess, 
but using LO prediction for OtG and NLO for SM leads to an 
artificial rise at large m t t- 


SM together with only LO prediction for O t G may lead to 
overestimate this effect, because the K factor decreases 
at larger energy scales. From Figure [5] we can see that 
the signal excess, a{O t G)/<- J (SM), is flat at a large energy 
range. In fact, at order 0( A -2 ), a(O t G) is suppressed 
by a constant factor of mf/A 2 compared with cr(SM), in¬ 
stead of s/ A 2 or y/smt/ A 2 as one might have expected 
naively by power counting. This is because the Higgs field 
in Otc always takes the vacuum expectation value, and 
the OtG operator flips the chirality of the top quark in 
its interference with the SM amplitude. Including higher 
order terms in 1/A 2 can give rise to additional contribu¬ 
tions that will indeed rise faster at large s, but if such 
an effect is large, it would imply the breakdown of the 
effective operator framework since the expansion in 1/A 
does not converge at large energy. In Table [Hl| we show K 
factors for SM and O t G as well as a(O t G)/ cr (SM ), with 
no cuts and with cuts m t i > 1 TeV and 2 TeV.Both LO 
and NLO give similar results for the signal excess, but 
using LO prediction for O t G with NLO for SM leads to 
an artificial rise at large m t i, due to the decreasing K 
factor of the SM. 


Another interesting observable is the forward- 
backward asymmetry, Apg, which has been observed at 
the Tevatron both by DO pS] and CDF [55]. Afb is de¬ 
fined as the asymmetry with respect to A y = y t — yi- In 
the SM the first non-zero contribution arises at NLO in 
QCD. In a dimension-six Lagrangian, only four-fermion 
operators can give a contribution at the tree level, and so 
Afb is another important observable that distinguishes 
between different new physics scenarios. In this respect, 
it is useful to know the first non-vanishing contribution 
from the CMDM operator, which appears at NLO, to at 
least have some estimation of the corresponding theoret¬ 
ical uncertainty related to this quantity. In our frame¬ 
work this calculation is straightforward. We expand the 
numerator and the denominator of Afb to NNLO for the 
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FIG. 5. Top quark pair invariant mass distribution above 1 
TeV, at the LHC 14 TeV energy. The renormalization and 
factorization scales are taken to be 1 TeV. The lower panel 
shows that the K factor for OtG decreases at high mass re¬ 
gion, and the signal over background ratio, cr(OtG)/a(SM), 
is almost a constant. 


SM part and NLO for the O t G part, and we find 

_ ^ew + Qig-^3 + oi^Ni + oigtj^NtG + tt s^ 2 ) 

FB a 2 s D 2 + a^D 3 + a*D 4 + aj^-D tG + 0(a 5 s ,a f A~ 2 ) 

=A fb (SM) + + 0(a 3 s ,a 2 s A~ 2 ) 

=0.095 ± 0.007 + C tG 0 . 021±°;™3 ^50 ( 14 ) 

where the SM prediction at NNLO in QCD is taken from 
Ref. [57], and the uncertainties of the second term come 
from scale variation. We thus expect a small modification 
to A F b from a non-vanishing CMDM. Given the current 
limit on CtOi however, this contribution is much smaller 
than the experiment uncertainties. 

Once differential cross sections are known, one can con¬ 
sider constraining the CMDM by using the normalized 
distributions of the tt observables. Unfortunately from 
Figure [3] and Figure [4] one can see that the shapes of 
the distributions from OtG are not significantly different 
than those from the SM. As a result the limits obtained 
only by using the shape of the distribution will be loose. 

As an example we consider the tt invariant mass dis¬ 
tribution at 7 TeV measured by the CMS collaboration 
[58] , We take only the first four bins, i.e. from 345 to 
650 GeV, to ensure the validity of the expansion in 1/A 2 
for A around TeV scale. We perform a simple x 2 fit f° r 
the differential cross section normalized within these four 
bins. We add the experimental and theoretical errors in 
quadrature. The SM prediction is computed with Mad- 
Graph5_aMC@NLO and then normalized to the most 
accurate NNLO+NNLL prediction |49| . with uncertain- 



FIG. 6. Transverse momentum distribution of the hardest 
muon at LHC 13 TeV. 


ties coming from the renormalization and factorization 
scale variation. The 95% CL allowed region is [-5.0,12.8], 
using LO predictions for the OtG contribution, and [- 
0.6,10.9] when using NLO prediction. Despite a signifi¬ 
cant improvement at NLO, the limit itself is much looser 
than those obtained from total cross sections, therefore 
we expect only a small improvement when the distribu¬ 
tion information is combined together with the total cross 
section in such analyses. 


B. Decayed top quarks 

The above results indicate that the kinematic observ¬ 
ables constructed from stable tt system are not very sen¬ 
sitive to the size of the top-quark CMDM. We therefore 
move on to include top-quark decays, where we expect 
that the decay products preserve the spin information of 
the top quarks, and thus can be more sensitive to the 
dipole structure in the operator. As an example we fo¬ 
cus on the dimuon channel, where both top quarks decay 
semileptonically into a 6 -quark and a muon. We use the 
MadSpin package [34] to decay the top quarks, so that 
the spin correlation at LO accuracy is preserved in the 
simulation. 

In our simulation we use the anti -k F algorithm for 
the jets with radius R = 0.5. The following cuts are 
imposed to mimic the environment of a real detector: 
Pt{ j ) > 30 GeV, \r](j)\ < 2.5, pt{£) > 20 GeV and 
rj(£) | <2.5, where j refers to jets and £ to muons. At least 
two jets, from which at least one containing a 6 -hadron, 
and exactly one pair of isolated muons are required. The 
isolation criteria is achieved by imposing a maximum 
value of 0.15 on the ratio of the scalar sum of p F of 
all hadronic tracks within A R = \J Aif + A (j) 2 < 0.3 
around the muon candidate, to the transverse momen- 
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FIG. 7. Transverse momentum distribution of the hardest 
6-jet at LHC 13 TeV. 



FIG. 8. Transverse momentum distribution of the second 
hardest 6-jet at LHC 13 TeV. 


This last distribution is particularly important for the 
measurement of the OtG operator, because it is sensi¬ 
tive to the spin correlation between the top quarks. The 
anomalous top-quark CMDM affects the spin correlation 
of the tt system m, and its effects have been searched 
for by the CMS collaboration f29j. The contribution from 
the O t Q operator to this distribution, expanded linearly 
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1.5 


S 23 o ig (nlo) 
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1.5 

IA<KU)I 
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FIG. 9. Difference in azimuthal angle between the two se¬ 
lected muons at LHC 13 TeV. Note that the K factor changes 
in a way that enhances the deviation of the OtG contribution 
from the SM. 


in CtG, takes the following form 

(1 da \ /1 da \ CtG / I da \ 

\ad\A(j)\) \a d\A<l)\ J SM + A 2 \cr d\A4>\) Np 


(15) 

which is valid provided C t o/3i/A 2 -C a$M- In the spirit 
of perturbation theory, the second term on the r.h.s. of 


Eq. (15) can be expanded to G(as)' 


1 da 


a d\A<p\ 


NLO 


NP 


1 da \ 
er d\A<f>\ ) 


LO 


NP 


1 da \ (1) 


NP 


turn of the muon. We show in Figure [6] the transverse 
momentum distribution of the muon. In Figure [7] [8] we 
show the hardest (largest pr) and the second hardest 
(when present in the event) 6-jets, respectively. In Fig¬ 
ure [9] the azimuthal angle difference between the two se¬ 
lected muons is shown. 


where 

/I da \ L ° 1 (do'gA 

\<t rf|A0|/ NP \d\^\ J 

and 


ffi° f d^SM \ 
^sm 2 \d|A</>|y 


(17) 


( 1 dcr V 1} = _JL ( da oto \ _L_ Lw ( da sM_\ 

\ad\A(t)\J Np a^yd\A<j)\) cr ^ 2 [ 1 \d\A<t>\) 


oho ( dc sm \ , (!) ( da o?G \ 

1 yd|A(/)| J SM yd|A^| J 


/ do-gM \ 

°sm 3 \d|A<^|y 
(18) 
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where dao tG represents the distribution from operator 
O t G f° r Ctc/A 2 = 1 TeV -2 . The superscript (1) indi¬ 
cates the as correction to the corresponding LO quan¬ 
tity. 


The predicted distribution is shown in Figure 10 for 
Ctc/A 2 = 1 TeV -2 . The OtG contribution has a pecu¬ 
liar structure which tends to flatten the distribution. In 


Figure [TT] we show the O t c distributions solely, as de¬ 
fined in Eqs. (16 1 (18). The first observation is that the 


purple and the red curves are very close to each other, 
indicating that the LO and NLO results are very simi¬ 
lar. The reason is that we are plotting the normalized 
distribution, and since the K factors are almost the same 
for the SM and for the O t G, they cancel each other when 
taking the ratio. In fact, one can see that Eq. van¬ 
ishes if the K factor is a constant. Alternatively, if we 
use NLO prediction for the SM but only LO prediction 
for O t G, following the logic that the radiative correction 
on the new physics effect is of higher order, then we will 
have 


<L 

!< 

;o 

O 


pp->tt at LHC13 
(N)LO+HERWIG6 

C tG /A 2 =1 TeV' 2 

— SM(NLO)+O tG (NLO) 

— SM(NLO)+O tc (LO) 

— SM(NLO) 

SM(LO) 


O 

z 


1.5 

N>(")l 


1 da 

a d| A(j)\ J 


\ nlo 


NP 


l 

t nlo 

7 sm 


( d °b? G \ 
\ d\A(j>\ J 


Pi° ( ^°SM \ 
a s N L02 \d\A<j,\) 

(19) 

This result, after expanding in a s , con tains only part of 
the 0(a s hr 2 ) corrections in Eq. (18) (and so we refer 
to as “nlo”). They come from the 0(a s ) corrections to 
the normalization, but not directly to the OtG contribu¬ 
tion. The missing 0(a s A -2 ) terms actually make a large 
difference, as illustrated by the blue curve in Figure El 
One can see that Eq. (19) gives a much lower estima¬ 
tion for the effect of OtG- This is not only because of 
the overall size of the K factor, but also due to the fact 
that the K factor is a decreasing function of A (f>(ll), and 
so the way it changes adds coherently to the difference 
in shapes between OtG and SM distributions, as can be 
seen in Figure [9j As a result, using NLO prediction for 
the SM together with only LO prediction for O t G signif¬ 
icantly underestimates the power of A <fi(ll) in discrimi¬ 
nating the O t G contribution from the SM. Also note that, 
the fact that Eq. (19) and the LO prediction in Eq. © 
differ implies that the LO prediction for OtG has a large 
uncertainty due to the missing 0(a s A -2 ) terms, which 
turn out to have a large effect in this special case. Thus 
our work improves the precision level of this prediction 
by completing the missing 0(a s Kr 2 ) terms. 

We also show for completeness two more angular dis¬ 
tributions which have been studied by Ref. m- Follow¬ 
ing Ref. m we define £~ (£ + ) as the momenta of the 
(anti-)muon in the rest frame of anti-top (top) quarks 

and k ( k ) the momenta of the top (anti-top) in the zero- 
mometum frame. In Figure [12] we show the distribu¬ 
tion of cos 9 1 cos 0 2 , where 9\ (# 2 ) is the angle Z(£~,k) 
(Z(£+,fc)). In Figure 


13 


we show the normalised distri¬ 
bution of cos#*, where 9* is the angle Z (£_,£+). Con¬ 
trary to the \A(f>(££)\ case, where the QCD corrections en¬ 
hances the anomalous coupling contribution to the shape, 


FIG. 10. Normalized distributions of the difference in az¬ 
imuthal angle between muons. 
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FIG. 11. New physics contribution to the normalized dimuon 
distribution at the LHC (13 TeV), da/d\A<j>\)^p defined 
in Eqs. fl5| fl9|). 


in these cases, we observe a uniform QCD correction with 
no effect in the normalised distributions. 
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FIG. 12. cos cos 02 distribution at LHC 13 TeV. 
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FIG. 13. cos0* distribution at LHC 13 TeV. 


VI. CONCLUSIONS 


The QCD correction increases the overall contribution 
from the OtG operator. For the total cross section for 
example, the increase is, at central scale, 12%, 43% and 
48% for Tevatron, LHC8 and LHC 14 respectively. More¬ 
over, the NLO calculation significantly reduces the scale 
uncertainty of the contribution from OtG • Limits on the 
coefficient of C t G are therefore improved. Our predicted 
allowed range at 95% CL using Tevatron and LHC8 data 
is —0.32 < CtG < 0-30 (assuming A = ITeV), which in 
terms of dy parameter gives —0.0096 < dy < 0.0090. 

Our implementation can be used for various exclusive 
studies. We have shown representative distributions for 
both stable and decayed top quarks as examples. We 
observed a significant reduction of scale variation in all 
distributions. The differential K factor is not a constant, 
but for all observables we have studied, it is similar to 
the SM K factor. Therefore we expect that using the SM 
K factor to rescale the LO contribution of OtG can be a 
good approximation for an NLO prediction in most cases. 
On the other hand, using NLO SM prediction together 
with LO prediction of OtG can be misleading in analysis 
where the ratio between O t c contribution and SM con¬ 
tribution can play a role. Observables sensitive to spin 
correlation can also be studied in the same framework, 
provided that the MadSpin package is used to preserve 
the spin information of the top quarks. This is particu¬ 
larly useful for spin correlation measurements where lim¬ 
its can be set by using various angular distributions of 
the decay products. We showed that the NLO correc¬ 
tion does not significantly change the LO prediction, but 
instead it increases the precision level, in particular for 
the A 4>(ll) distribution, where using NLO SM prediction 
together with LO prediction of O t G can lead to underes¬ 
timate the effect from top-quark CMDM. 

Our theoretical approach is based on the effective field 
theory for top-quark couplings, and is a first step of the 
automation of the top-quark flavor-diagonal operators 
in the MadGraph5_aMC@NLO framework. The next 
step is to extend our study to other top-quark operators, 
including the CP-odd ones such as the CEDM, as well as 
other electroweak couplings of the top quark. These stud¬ 
ies will pave the way to a global analysis for top quark 
couplings using the effective field theory framework. 


In this work we have presented the NLO calculation 
for top-quark pair production, including an anomalous 
top-quark CMDM, as described by the dimension-six op¬ 
erator OtG ■ Our calculation is implemented in the Mad- 
Graph5_aMC@NLO framework, which allows the re¬ 
sult to be matched to parton shower automatically. We 
have studied the impact of QCD corrections to the contri¬ 
bution of the CMDM in top quark pair production, both 
on total cross section as well as on various distributions. 
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